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j>! . ABSTRACT 

\^ ' We present four X-ray/7-ray spectra of Cyg X-l observed in the hard ('low') state 

I/"") I simultaneously by Ginga and GRO OSSE on 1991 June 6. The four spectra have 

almost identical spectral form but vary in the normalisation within a factor of two. 
The 3-30 keV Ginga spectra are well represented by power laws with an energy spectral 
index of a ~ 0.6 and a Compton reflection component including a fluorescent Fe Ka 
corresponding to the solid angle of the reflector of ~ 0.3 x 2n. These spectra join 
smoothly on to the OSSE range (> 50 keV) and are then cut off above ~ 150 keV. 
. The overall spectra can be modelled by repeated Compton scattering in a mildly- 

1' relativistic, thermal plasma with the optical depth of r ~ 1. However, the high-energy 

JL \ cutoff is steeper than that due to single-temperature thermal Comptonisation. It can 

^**| ■ be described by a superposition of dominant r ~ 1-2, thermal emission at kT ~ 100 

keV and a Wien-like component from an optically-thick plasma at kT ~ 50 keV. 
, The X-ray spectra do not show the presence of an anisotropy break required if 

thermal Compton scattering takes place in a corona above a cold disc. Also, the flat 
spectral index shows that the plasma is soft-photon starved, i.e., the luminosity in 
, incident soft X-ray seed photons is very much less than that in the hard X-rays. Fur- 

5_i ■ thermore, the observed solid angle of the reflector is significantly less than 2tt. These 

facts taken together strongly rule out a disc-corona geometry. Rather, the observed 
spectra are consistent with a geometry in which the cold accretion disc (which both 
supplies the seed soft X-rays and reflects hard X-rays) only exists at large radii, while 
the Comptonising hot plasma is located in an inner region with no cold disc. This hot 
plasma consists of either pure e ± pairs if the source size is ~ 5 Schwarzschild radii or 
it contains also protons if the size is larger. 

Key words: accretion, accretion discs - gamma-rays: observations - gamma-rays: 
theory - stars: individual (Cygnus X-l) - X-rays: stars 



1 INTRODUCTION 

Cyg X-l, one of the brightest X-ray sources in the sky, was 
discovered in X-rays more than 30 years ago (Bowyer et al. 
1965). This binary system is at a distance of ~ 2.5 kpc and 
consists of the 09.7 lab type supergiant HDE 226868 (Gies 
& Bolton 1986) orbiting around a compact object with a 
period of 5.6 days. The mass of the unseen companion is 
significantly larger then 5 Mq (Dolan 1992), strongly sug- 



gesting that it is a black hole. Focused wind accretion (Gies 
& Bolton 1986) from the primary star (which is extremely 
close to filling its Roche lobe) drives the powerful source 
of the X-ray radiation. The X-ray emission of Cyg X-l ex- 
hibits strong variability on all time scales from milliseconds 
to months (e.g. Miyamoto & Kitamoto 1989; Priedhorsky, 
Terrel & Holt 1983). However, it spends most of the time 
in the hard (so-called 'low') X-ray state, characterised by a 
low flux in soft X-rays and strong hard X-ray and soft 7-ray 
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fluxes (e.g., Phlips et al. 1996, hereafter P96). In this state, 
the 50 keV flux changes by about an order of magnitude but 
the spectral variability is modest (P96). From time to time, 
Cyg X-l is in the soft ('high') state (e.g., Liang & Nolan 
1984, Cui et al. 1997), with a strong soft X-ray emission 
and a relatively weak tail in hard X-rays. 

The hard-state spectrum of Cyg X-l at i; 3 keV con- 
sists typically of a power law with an energy spectral index 
of a ~ 0.6-0.7 and a Compton-reflection continuum compo- 
nent, prominent above 10 keV (Done et al. 1992; Ebisawa et 
al. 1996, hereafter E96). This component also includes an Fe 
K edge above ~ 7 keV and an Fe Ka fluorescence line (Barr, 
White & Page 1985; Marshall et al. 1993; E96). The results 
of BBXRT (Marshall et al. 1993) and ASCA (E96) indicate 
that the line is narrow, with the width of a ;S 0.2 keV, sug- 
gesting its origin in outer parts of an accretion disc. The 
edge energy suggests the reflecting surface is mildly ionised, 
with Fe XI as the most abundant Fe species (E96) . 

The soft 7-ray spectra of Cyg X-l in the hard state are 
steeply cut off above ~ 150 keV (P96). P96 obtained good 
fits to the GRO OSSE spectra of Cyg X-l with exponen- 
tially cut off power laws. However, the X-ray spectral index 
implied by those fits is much harder than that observed in 
X-rays, showing that the broad-band X-ray /7-ray (hereafter 
X7) spectrum is more complex. 

The most likely physical process responsible for the in- 
trinsic continuum appears to be thermal Comptonisation of 
soft X-rays from thermal emission of a cold accretion disc at 
a temperature of ~ 100 eV (e.g., Balucinska-Church et al. 
1995). Haardt et al. (1993) pointed out that Comptonisation 
may take place in an optically thin plasma rather than in an 
optically thick one (as fitted by, e.g., Sunyaev & Triimper 
1979). 

In this work, we present four 2-1000 keV spectra of Cyg 
X-l observed in the hard state simultaneously by Ginga and 
OSSE on 1991 June 6. This provides us with an unprece- 
dented opportunity to study physical processes taking place 
in the X7 source as well as the source geometry. We consider 
thermal and nonthermal Comptonisation, Compton reflec- 
tion, e pair production, and constraints on the sizes and 
relative covering factors of the X7 source and cold matter. 



2 THE DATA 

Cyg X-l was observed four times by Ginga on 1991 June 6 
during the OSSE viewing period 2 (P96). Overlapping ob- 
servation periods provide us with four sets of nearly simulta- 
neous data. The orbital phase during the four observations 
was from 0.70 to 0.85. Table 1 gives the log of observations. 
We adopt 4661 cm 2 as the Ginga effective area, which is the 
standard value used in the Leicester data base (D. Smith, 
private communication) . The high voltage of the Ginga LAC 
detector (Turner et al. 1989) was reduced in this observa- 
tion, which lead to a usable energy range of 2-30 keV (E96) . 
A 1 per cent systematic error was included in each Ginga 
channel. 

The OSSE data are from 50 keV to 1000 keV. They 
include estimated systematic errors computed from the un- 
certainties in the low energy calibration and response of the 
detectors using both in-orbit and prelaunch calibration data. 
The energy-dependent systematic errors are expressed as an 
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Figure 1. The four spectra of Cyg X-l, labelled by the obser- 
vation number (see Table 1), observed simultaneously by Ginga 
and OSSE on 1991 June 6. The solid curves represent the two- 
temperature model described in Section 3.1. The data have been 
rebinned for clarity of the display. 

uncertainty in the effective area in the OSSE response. These 
systematic errors were added in quadrature to the statistical 
errors. The former are most important at the lowest energies 
(approximately 3% uncertainty in effective area at 50 keV, 
decreasing to 0.3% at 150 keV and above). 

Figure 1 shows the four unfolded spectra corresponding 
to the four data sets. (The solid curves represent a two- 
component model described in Section 3.1 below.) The ra- 
tio of the highest 2-1000 keV model flux to the lowest one, 
obtained during the second and the fourth observation, re- 
spectively (separated by about 15 hours), is ~ 1.8. 



3 RESULTS 
3.1 Spectral fits 

In our fits, we use the xspec spectral fitting package version 
9 (Arnaud 1996). Our model spectra are absorbed by a col- 
umn density, Ah, which we constrain to be > the Galactic 
column, 6±2 x 10 21 cm~ 2 (Balucihska-Church et al. 1995). 
We obtain absorption larger than the Galactic one for some 
models, especially for the data sets 3 and 4 (Tables 2 and 3), 
which can be attributed to transient dips associated with the 
intervening gas (Kitamoto et al. 1984). In X-rays, the spec- 
tra consist clearly of two components: an underlying power 
law and a component due to reflection from the surface of an 
accretion disc including both continuum Compton reflection 
and a fluorescent Fe Ka line (e.g., George & Fabian 1991). 

We first fit the Ginga data only with the incident con- 
tinuum in the form of a power law e-folded with E c = 300 
keV together with the Compton reflection component and 
an iron line. In the spectral fitting we ignore the data below 
3 keV due to the presence of a soft excess at lower energies 
(E96), but we show these data in the plots (Figs. 1-4). We 
use angle-dependent Compton reflection Green's functions 
of Magdziarz & Zdziarski (1995). The disc inclination is as- 
sumed to be less than 67° (Dolan & Tapia 1989), and in 
most fits the most probably inclination of 30° is used (Gies 
& Bolton 1986). The luminosity intercepted by the reflecting 
medium equals Q,/2-k times the luminosity emitted outward 
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Table 1. The log of observations in UT on 1991 June 6. The counts of Ginga and OSSE arc for 2-30 keV and 50-150 keV, respectively. 





Ginga 








OSSE 








Data set 


Start 


End 


Live time [s] 


Counts/s 


Start 


End 


Live time [s] 


Counts/s 


1 


00 h 17 m 50 s 


02 h 09 m 50 s 


2304 


3410 ± 10 


00 ,, 02 m 32 s 


00 ,, 37 m 29 s 


1806 


455 ±4 












01 h 23 m 21 s 


02 h ll m 09 s 


2315 


554 ±4 


2 


04 h 43 m 26 s 


06 h 29 m 02 s 


888 


3890 ± 10 


04' l 28 m 46 s 


05' l 03 m 43 s 


1741 


456 ±5 












06 h 02 m 26 s 


06 h 50 m 29 a 


2299 


665 ±5 


3 


ll h 02 m 38 s 


14 h 24 m 46 s 


2828 


2290 ± 10 


10 h 42 m 50 s 


11 h 15 m 36 s 


1708 


400 ±4 












12 ft 14 m 47 s 


12 h 54 m 07 s 


2010 


371 ±4 












13 h 48 m 41 s 


14 h 32 m 23 s 


2257 


307 ±4 


4 


20 h 22 m 14 s 


20 h 43 m 58 s 


1272 


2080 ± 10 


20 h 02 m 15 s 


20 ,l 32 m 50 3 


1629 


316 ±4 



by the X7 source, where Q is the solid angle subtended by 
the reflector as seen from the source of isotropic radiation. In 
some models, the reflecting medium is allowed to be ionised 
with an ionisation parameter, £ = L/(nr 2 ), where L is the 5 
eV-20 keV luminosity in a power law spectrum and n is the 
density of the reflector located at distance r from the illumi- 
nating source (Done et al. 1992). The reflector temperature 
is kept at 130 eV, which corresponds to the temperature of 
the soft X-ray spectral component (Baluciriska-Church et al. 
1995). The abundances are from Anders & Ebihara (1982) 
except that the relative Fe abundance, Af c , is a free param- 
eter. The ion edge energies and opacities are from Reilman 
& Manson (1979) except for the Fe K-edge energies taken 
from Kaastra & Mewe (1993). 

We find little intrinsic spectral variability in the Ginga 
data, so we use the average Ginga spectrum to better con- 
strain the parameters in X-rays. We obtain a — 0.59l°'°,3, 
U/2-K = 0.34l°;°|?, A Fc = 2.1±g;| ( X 2 = 21/22 d.o.f.). The 
ionisation of the reflector is weak, £ = 23^22 er 8 cm s_1 - 
There is a distinct Fe Ka line in the spectrum, which we 
model here as a Gaussian centred at 6.4 keV and with the 
width of 0.1 keV, consistent with both the Ginga resolution 
not allowing us to resolve the line and the ASCA results in- 
dicating the line is narrow (E96) . The line equivalent width 
is EW=74^22 eV. These results are similar to those obtained 
for other Ginga observations (e.g. E96). 

In the fits with both Ginga and OSSE below, we exam- 
ine the effect of allowing the relative normalisation of the 
two data sets to be free within ±15 per cent (which is a con- 
servative limit to the relative calibration uncertainty of the 
two instruments). However, we find the relative normalisa- 
tion to be consistent with unity as allowing it to be free leads 
to marginal fit improvements only. Therefore, in tables and 
figures below we use the actual normalisation of the Ginga 
and OSSE spectra. 

When we extrapolate the Ginga model to the OSSE 
range, there is a good match around 50 keV, but the OSSE 
data require a cutoff in the spectrum at J; 150 keV. How- 
ever, a simple e-folded power law together with its reflection 
spectrum gives a cutoff much too shallow to account for the 
sharp break seen in the data, e.g., xt = 220/76 for the data 
set 2. (The best fit parameters for this data set are: a — 0.52, 
E c = 270 keV, Q/2ty = 0.28.) 

However, e-folded power laws approximate only roughly 
thermal Comptonisation spectra, as discussed by, e.g., 
Poutanen & Svensson (1996; hereafter PS96). Therefore, we 
then use a thermal-Compton disc-corona model of PS96. 
In the model used, the hot coronal plasma is assumed to 



form hemispheres (with the radial Thomson optical depth 
of r) on the surface of an accretion disc. This model, how- 
ever, is also very strongly ruled out by the data, yielding 
xt ~ 1000/75. The main reason for the bad fit is the pres- 
ence of an anisotropy break in the underlying model spec- 
trum. The anisotropy break appears in the disc-corona spec- 
trum viewed from above (at i = 30°) around the peak en- 
ergy of the second-order scattering because the first-order 
scattering is directed mainly towards the disc (Stern et al. 
1995; PS96; Svensson 1996). For the seed blackbody pho- 
tons at the temperature of 130 eV (Balucinska-Church et 
al. 1995) and the hot plasma temperature of ~ 100 keV, the 
anisotropy break energy is ~ 6 keV. The Ginga data are 
instead very well described by a single power law up to 30 
keV and a Compton-refiection component on top of it, as 
found above. A secondary problem with the model of PS96 
applied to the present data is the predicted solid angle of the 
reflector of Q = 27r, whereas the data prefer a lower value. 
Still, allowing £1 < 2n in the model (which may correspond 
to a cold disc present only at a limited range of radii) does 
not result in an acceptable fit. E.g., the fit to the data set 
1 (Fig. 2) yields %l = 160/75 for kT = 92 keV, r = 2.2, 
0/27T = 0.24 for i = 30°. The model can be strongly ruled 
out (x 2 = 117) even for the largest viewing angle allowed by 
the optical data, i = 67° (Dolan & Tapia 1989) The reason 
for the bad fit of the corona model is still the presence of 
the anisotropy break in the incident model spectrum. E.g., 
the incident spectrum in Figure 2 (dashed curve) can be ap- 
proximated by a harder power law, a ~ 0.43, below 5.8 keV 
and a softer one, a ~ 0.56 above it. The pattern of residu- 
als in Figure 2 clearly show how this broken power-law-like 
model does not fit the data (which are consistent with a 
single incident power law). 

We have also considered the homogeneous disc-corona 
model of Haardt & Maraschi (1993) as a phenomenological 
description of the data (although that model does not satisfy 
the energy balance in the case of Cyg X-l, see Section 3.2). 
For that model, we use the code of Haardt (1993). The main 
difference with respect to the hemisphere model above is 
that now there is strong Comptonisation of the reflection 
component by the hot corona with r ~ 1 above the disc 
(see, e.g., PS96). Again, no acceptable fit is obtained, xl ~ 
200/75. Thus, the spectral data of Cyg X-l rule out the 
disc-corona model regardless of the corona geometry. 

Motivated by our ruling out the disc-corona model, we 
have considered thermal Comptonisation in a hot plasma 
cloud irradiated approximately isotropically by soft seed 
photons. This removes the presence of the anisotropy break 
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Figure 2. The anisotropic disc-corona model (PS96) fitted to the 
data set 1. The upper panel shows the data (crosses), the ther- 
mal Comptonisation continuum (dashed curve) , and the reflection 
component (dotted curve). The solid curve shows the sum. The 
lower panel in this figure and in figures below shows the contri- 
bution to the total x 2 from separate data channels [multiplied by 
the sign of (data — model)]. This model provides no satisfactory 
fit to the data due to the incident spectrum having a break at 
~ 5 keV. 



in the spectrum, but it retains the proper description of 
the high-energy cutoff from Comptonisation. To this end, 
we used an isotropic-scattering version of the Comptonisa- 
tion code of PS96 (J. Poutanen, private communication). 
We have assumed that a fraction of 0.3 of the reflection 
spectrum undergoes further thermal Comptonisation in the 
hot cloud, which is consistent with source geometry inferred 
from energy balance (Section 3.2 below); the fit results are 
insensitive to this number as long as it is <C 1. Also, the fits 
are insensitive to the viewing angle; if i is increased from 
our default value of 30° to 67° , \ 2 remains the same within 
±1 for the four data sets. This model does not allow the re- 
flector to be ionised, which may result in some overestimate 
of x 2 since the Ginga data indicate the reflector is weakly 
ionised. The Fe abundance of this model consistent with all 
4 data sets is Af c = 1. As expected, we obtain a dramatic 
improvement of the fits, with \ 2 1 now. The fit results 
for the four data sets are given in Table 2, and the fit to 
the data set 2 is presented in Figure 3. We see that the hot 
plasma parameters are kT ~ 100 keV and r ~ 2 for all data 
sets. 

We have tested the isotropic Comptonisation code of 
PS96 by our own Monte Carlo calculations (see Zdziarski, 
Johnson & Magdziarz 1996). We assumed a spherical source 
illuminated by seed photons at the center. We obtain the 
same plasma temperature, kT ~ 100 keV, but the radial 
optical depth of the sphere is r ~ 1.3. The difference be- 
tween this t and r ~ 2 obtained above is explained by the 
difference between the spherical source in the Monte Carlo 



Figure 3. The isotropic one-temperature thermal Comptonisa- 
tion model fitted to the data set 2. The upper panel shows the 
data (crosses), the thermal Comptonisation continuum (dashed 
curve) and the reflection component (dotted curve). The solid 
curve shows the sum. This model yields xi < lj °ut there is 
still a systematic residual pattern for the OSSE data (the lower 
panel). 



calculations and the hemisphere source in the method of 
PS96. 

Although xt < 1 for the isotropic Comptonisation 
model, there are still strong systematic residuals in the fit to 
the OSSE data, see the lower panel of Figure 3. Specifically, 
the model spectrum cuts off at high energies less than the 
data, i.e., the model spectrum is broader than the observed 
spectrum. Such an effect could appear if there were an ad- 
ditional narrow spectral component at high energies. Such 
a sharp spectral component can be due to a hot, optically- 
thick plasma component in the source, e.g., a transition re- 
gion between the hot and the cold medium in the accretion 
flow (Misra & Melia 1996). 

We have therefore considered models with two hot- 
plasma components. The first component, along with its re- 
flection, is computed in the same way as in the isotropic 
Comptonisation model above. The Fe abundance is kept at 
Af = 1.0, which is both the average value from fits to the 
four spectra and a value within the confidence limits of all 
of the four fits. The second, additional component (whose 
reflection is not considered) is due to optically-thick thermal 
Comptonisation and it is computed using a modified Kom- 
paneets equation as in Zdziarski et al. (1996). We see in Ta- 
ble 3 and Figure 4 that addition of a second Comptonisation 
component with kT' ~ 50 keV improves significantly the fit. 
In the actual fits, we have had to prevent the secondary com- 
ponent from giving a dominant contribution at low energies. 
This is achieved by fixing its low-energy spectral index to 
a' = 0.15, which we found to have a negligible effect on \ 2 ■ 
The Thomson optical depth of the second plasma component 
has been found by the Monte Carlo method. It is r' ~ 6 for 
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Table 2. The parameters of the fits with the single-temperature isotropic model. / is the 1-keV normalisation in cm~ 2 s -1 , kT is in 
keV, and EW is in eV. The column density, JVjj (in units of 10 22 cm -2 ), is constrained to be > the lower limit on the Galactic column. 
The errors are given for 90 per cent confidence intervals, A\ 2 = 2.7 (Lampton et al. 1976), and \ 2 ls given for 76 d.o.f. 



Obs. 


N H 


kT 


T 


Q/2tt 




EW 


x 2 


1 


0.4+0- 14 


103± 7 5 


1 98+ 09 


34+0 03 


1.19 


1J ' -29 


69 


2 


04 +o.09 






32+ - 04 
u -° -0.03 


1.34 


137+ 29 


97 


3 


1 iq+0.18 


ioii 8 8 


2.08±°;» 


19+ 003 

u -0.03 


0.79 


121+ 3 " 


60 


4 


o.78t°,; 2 2 ? 


120t 2 g 


1.741°;- 


25+ - 04 


0.74 




64 



Table 3. The parameters of the fits with the two-component thermal-Comptonisation model. The parameters of the main component 
are denoted by the same symbols as in Tabic 1 , and I' and kT' are the 1-keV normalisation and temperature of the secondary component 
in units of cm~ 2 s _1 and keV, respectively. L is the 2-1000 keV luminosity in 10 37 erg s _1 for the assumed distance of 2.5 kpc, and \ 2 
is given for 74 d.o.f. 



Obs. 




kT 


T 


U/2tt 




kT' 


/' 


EW 


L 


x 2 


1 

2 
3 


0.85+ ,; 2 ? 
0.95+°; 27 
1 50+ ' 29 


104+1 7 
107+1? 
108+1 4 . 


1 76+0- 32 

1 - ' -0.26 

1 67+ ' 29 

1 -°'-0.28 


49+ 07 
"«+ - 07 
25+ 06 


1.33 
1.53 
0.83 


51+.?, 

47 +6 

43+21 


0.0080 
0.0106 
0.0026 


iik+31 
li0 -30 
1 1 0+ 32 

104+ 32 


3.7 
4.1 
2.6 


48 
64 
56 


4 


1 18+ 028 


127+ 27 


1.53lg;ig 


32+ 08 

u -° -0.06 


0.79 


36+ 13 

J -10 


0.0042 


S7 -ll 


2.3 


54 



a sphere with a central source of seed photons. The lumi- 
nosity in the best-fit optically-thick component is 13, 15, 5, 
and 8 per cent of the total luminosity for the data sets 1, 
2, 3, and 4, respectively. This model provides excellent fits 
to the data (\ 2 ~ 50-60/74 d.o.f.), and the improvement 
with respect to the one-temperature model is statistically 
significant at more than 99.9 per cent (the probability that 
the reduction of \ 2 is by chance is ~ 10~ 13 ) taking into ac- 
count all four data sets. The two-component model strongly 
favours a face-on geometry; when i is increased from 30° 
(assumed above) to i = 67° , A\ 2 — +7 for each of the data 
set 1 and 2. 

A possible alternative explanation to the hot plasma 
with a distribution of kT and r is Comptonisation by a 
plasma with nonthermal electron distribution. We have per- 
formed a number of Monte Carlo computations using ei- 
ther a Maxwellian electron distribution truncated at high 
energies or electron power-law distributions constrained to 
low values of the Lorentz factor 2). Our conclusion is 
that those models can in principle provide high-energy cut- 
offs stronger than those from the Maxwellian. However, \ 2 
fitting of Monte Carlo results would require a very exten- 
sive grid of models. Instead, an iterative-scattering method 
(PS96) modified for Comptonisation by a non-Maxwellian 
distribution appears to provide a simpler solution to this 
problem (J. Poutanen, private communication). We leave 
detailed investigation of this issue for future work. 

3.2 Energy balance, geometry, and pair 
production 

In Section 3.1 above, we have ruled out the disc-corona 
model for Cyg X-l based on spectral fitting alone. This con- 
clusion can be further supported by considering the energy 
balance between the hot and the cold phase. 

The presence of Compton reflection shows there is cold 
matter in the vicinity of the X7 source. Only some of the in- 
cident X7 photons get reflected (with the integrated albedo 
of ~ 0.15, e.g., Magdziarz & Zdziarski 1995), but most 
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Figure 4. The two-temperature thermal Comptonisation fit to 
the data set 1. The dashed and dot-dashed curves represent the 
primary and secondary Comptonisation components, respectively. 
The dotted curve represents Compton reflection of the primary 
component. The solid curve gives the sum. This model provides 
a statistically satisfactory description of the data. 



(~ 0.85) are absorbed and reemitted in soft X-rays with 
a distribution close to a blackbody. Depending on the ge- 
ometry, a fraction of the blackbody photons will return to 
the X7 source and provide some of the seeds for thermal 
Compton upscattering, which cools the hot plasma. (Fur- 
ther cooling can be provided by a blackbody flux from dissi- 
pation of gravitational energy in the cold matter). The ratio 
of the luminosity in the X7 spectrum to that of the seed 
photons (the Compton amplification factor) determines the 



© 0000 RAS, MNRAS 000, 000-000 



6 M. Gierlinski et al. 



resultant spectral shape because, at a given optical depth, 
the electron temperature is determined by the amount of 
cooling photons incident on the hot plasma. From Monte- 
Carlo simulations, we find that the rather flat continuum 
spectrum (a ~ 0.6) extending out to ~ 150 keV can only be 
produced if the incident seed photon luminosity is about 14 
times less than the luminosity in the Comptonised X7 spec- 
trum. This rules out a homogeneous disc-corona geometry 
of the source, in which approximately half of the hard X-ray 
photons are intercepted by the disc, 85 per cent of which are 
thermalised, and then re-intercepted by the corona, giving 
an X-ray luminosity only ~ 2.4 times that of the seed pho- 
tons. The lack of soft photons could in principle be obtained 
if the hot corona is in the form of detached active regions 
(e.g., PS96). However, this would imply an observed solid 
angle covered by the reflecting medium as seen from the hot 
source close to 2tt whereas the observed value is ~ 0.4 x 2tv. 
This means that the hot source does not form a corona (ei- 
ther homogeneous or patchy) above the surface of a cold 
disc, fully consistent with the results of spectral fitting in 
Section 3.1. 

Assuming that the reflection spectrum is not strongly 
Comptonised, i.e., the observed amount of reflection is in- 
dicative of the solid angle subtended by the cold material 
to the hard X-ray source, then we can estimate the frac- 
tion, g, of the reprocessed radiation that returns to the hot 
source to form the Compton seed photons. From energy 
balance considerations g < (0.85 x 0.5 x Q/2n x 14) _1 ~ 
(6 x D./2TT)- 1 ~ 0.4. (A similar value was used in the fits in 
Section 3.1.) Again this highlights the problems of the homo- 
geneous corona-disc geometry, in which all the reprocessed 
photons are re-intercepted by the hot plasma. 

The considerations above explain the fit results in which 
we have been unable to obtain satisfactory models with the 
disc-corona model (see above). A possible source geometry 
consistent with the observed spectra is a hot inner region 
together with a colder outer disc (e.g. Shapiro, Lightman & 
Eardley 1976; Bjornsson et al. 1996; Narayan 1996). This 
geometry explains the observed low reflection fraction, the 
large Compton amplification of seed photons, as well as the 
narrowness of the Ka line observed by E96. Also, the seed 
soft photons in this model are not produced underneath the 
hot plasma but rather come from various directions side- 
ways. This then explains the absence of the anisotropy break 
in the spectrum. We intend to further investigate the ex- 
pected spectrum from such geometry in future work. 

The last issue we consider is the presence of pairs in 
the hot source. Knowing the shape of the spectrum and r 
(~ 1-2), we can compute the compactness parameter, I = 
Lot / Rrn c c 3 (where R is the source size), at which photon- 
photon pair production rate is sufficient to produce pairs 
with the given r (Svensson 1984). We obtain i ~ 40-70 
roughly correlated with the luminosity varying in the range 
of ~ 0.02-0.03 of the Eddington limit (assuming the mass of 
10M©). This implies an approximately constant lower limit 
on the source size of R ~ 5 Schwarzschild radii. The source 
is pair-dominated at this limit, and is electron-dominated 
if R is much larger. The constancy with L is due to an 
approximate constancy of the model fluxes seen around 0.5 
MeV (in spite of the overall X7 flux varying), see Figure 1. 

We also note that the pair annihilation feature at the kT 
and r implied by the thermal Comptonisation model with 



e pairs would be rather weak, with the peak annihilation 
flux less than 0.1 of the observed 0.5 MeV flux (see Maciolek- 
Niedzwiecki, Zdziarski & Coppi 1995), and the luminosity in 
annihilation photons of less than ~ 2 per cent of the total 
luminosity. Thus, the observed spectra are consistent with 
(but do not require) the hot (kT ~ 100 keV) source made 
of e pairs. 

Note that virtually all pair production is due to the 
hotter and optically thinner plasma component, with kT ~ 
100 keV and r ~ 1-2. This can be seen in Figure 4, in 
which number of Wien photons around 0.5 MeV is much less 
than that in the main spectral component. Consequently, 
the cooler, optically thick, plasma (kT ~ 50 keV, r ~ 6) is 
made entirely of electrons and protons. 



4 CONCLUSIONS 

Wc find the X7 source in Cyg X-l does not form a corona 
above the surface of an accretion disc. Three arguments lead 
us to that conclusion: the solid angle subtended by the cold 
medium as seen from the X7 source is significantly smaller 
than 2-k, the X7 source is photon-starved, and there is no 
anisotropy break (necessary in the disc-corona models) in 
the X-ray spectrum. In addition, the reflecting cold medium 
is far away from the central black hole, as demonstrated by 
the narrowness of the fluorescent Ka line observed by ASCA 
(E96). All the properties above are consistent with a hot 
inner disc Comptonising soft photons from a cold outer disc. 
This distinguishes Cyg X-l from some Seyfert Is, in which 
the cold, reflecting, disc appears to extend all the way to the 
innermost stable orbit (e.g., Fabian et al. 1995). Note that 
this difference cannot be explained by the radiation-pressure 
disc instability in a disc-corona system, which is stronger for 
more massive black holes (e.g., Svensson & Zdziarski 1994). 

The most likely radiative process producing X7 pho- 
tons in Cyg X-l is Comptonisation of seed soft X-rays 
photons. The observed sharpness of the high-energy cutoff 
rules out a uniform thermal plasma. Either a distribution 
of plasma parameters or an electron distribution different 
from a Maxwellian is required. A remarkable property of 
the spectra is their almost constant shape with the vary- 
ing amplitude. In models with repeated Compton scatter- 
ing, this implies the constancy of the plasma parameters 
(e.g., t and kT for thermal Comptonisation) with changing 
accretion rate. This represent a strong constraint on future 
models of Cyg X-l. In particular, simplest thermal e ± pair 
models predict a strong dependence of r on the luminosity 
(e.g. Zdziarski 1985). The constant spectral shape can be 
achieved in models with high-energy nonthermal electrons 
(e.g., Svensson 1987; Lightman & Zdziarski 1987). However, 
these models predict a strong annihilation feature which is 
not seen in the data (P96), although this constraint could 
be possibly weakened if the annihilation rate were reduced 
by having the pairs reaccelerated before they lose energy 
completely. 
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